A Well-balanced Finite Volume-Augmented Lagrangian Method for an Integrated Herschel-Bulkley Model

We are interested in the derivation of an integrated Herschel-Bulkley model for shallow flows, as well as in the design of a numerical algorithm to solve the resulting equations. The goal is to simulate the evolution of thin sheet of viscoplastic materials on inclined planes and, in particular, to be able to compute the evolution from dynamic to stationary states. The model involves a variational inequality and it is valid from null to moderate slopes. The proposed numerical scheme is well balanced and involves a coupling between a duality technique (to treat plasticity), a fixed point method (to handle the power law) and a finite volume discretization. Several numerical tests are done, including a comparison with an analytic solution, to confirm the well balanced property and the ability to cope with the various rheological regimes associated with the Herschel-Bulkley constitutive law.